Shot noise in a diffusive F-N-F spin valve 
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Fluctuations of electric current in a spin valve consisting of a diffusive conductor connected to 
ferromagnetic leads and operated in the giant magnetoresistance regime are studied. It is shown 
that a new source of fluctuations due to spin-flip scattering enhances strongly shot noise up to a 
point where the Fano factor approaches the full Poissonian value. 
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Transport in various spintronic devices contain- 
ing ferromagnet-paramagnet interfaces is attracting a lot 
of attention. Considerable experimental and theoreti- 
cal efforts have been directed towards the understanding 
of magnetoresistance, spin injection, spin accumulation, 
spin-orbit interaction, current-induced torque and other 
fascinating and challenging effects (the vast and quickly 
expanding bibliography is far beyond the scope of this 
Letter). Advances in technology and sample fabrication 
resulting in devices of nanoscale dimensions led the meth- 
ods and notions of spintronics to be the natural outgrows 
and further developments of the exciting and successful 
ideas of mesoscopics. 

One of the issues outstanding in mesoscopic physics has 
been the phenomenon of the shot noise,i.e. current fluc- 
tuations in nonequilibrium conductors |2f. In particular, 
an experimental confirmation [3j of the theoretically pre- 
dicted 1/3-suppression (compared to the Poissonan value 
characteristic for the transmission of independent parti- 
cles) of the noise signal in diffusive conductors is one 
of the milestones in the field. Shot noise in ferromagnet- 
normal metal constrictions is also evolving into a subject 
of much interest. Current fluctuations in a F-quantum 
dot-F system in the Coulomb blockade regime were con- 
sidered in Refs. 0, II 110, 

noise in a quantum dot in 
the Kondo regime analyzed in Ref. 10], ballistic beam 
splitter with spin-orbit interaction discussed in Ref. . 
Dependence of the shot noise in a diffusive conductor at- 
tached to ferromagnetic reservoirs on the relative angle 
between the magnetizations of reservoirs has been stud- 
ied in Ref. 01 with the help of the circuit theory [l3j| . 
However, effects of a spin-flip scattering on the fluctua- 
tions of electric current in diffusive conductors have been 
disregarded so far. In the present Letter we show them 
to make a profound effect on the shot noise power. 

The universal 1 /3-shot noise in a conventional diffusive 
conductor is due to the interplay of the random impurity 
scattering and restrictions imposed by the Fermi statis- 
tics. In the presence of ferromagnetic contacts, however, 
the spin degeneracy is lifted with spin-up and spin-down 
electrons representing two different subsystems. The 
number of particles in each subsystem is not conserved 
(due to spin-flip scattering) leading therefore to a new 
class of fluctuations. The situation here resembles closely 



the fluctuations of radiation in random optical media |14j . 
The absence of particle conservation in a gas of photons 
results in the enhancement of photon flux noise above the 
Poissonian value (also the result of bunching typical for 
bosons). With the notable difference in statistics (Fermi 
instead of Bose) the framework of stochastic diffusion 
equations [lil Hq can be formulated for the fluctuations 
in disordered spintronic devices as well. 
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FIG. 1: Spin valve consisting of a paramagnetic diffusive con- 
ductor (N) connected to ferromagnetic leads (F) through tun- 
nel contacts. Conduction bands in the leads are assumed to 
be completely spin-polarized. For small amount of spin-flip 
scattering the resistance to electric current is large when the 
magnetizations are antiparallel ('off '-state of the valve), com- 
pared to the usual metallic resistance for the parallel config- 
uration ('on '-state). 

To demonstrate this we discuss the most characteristic 
example of a spin valve in the giant magnetoresistance 
regime, when the transport across the valve is extremely 
sensitive to the intensity of a spin-flip scattering. Namely, 
we consider a diffusive paramagnetic conductor (N) sand- 
wiched between two ideal ferromagnetic (F) leads, Fig. 1. 
'Ideal' means that electron distributions inside the leads 
are not affected by the presence of the normal region (a 
typical mesoscopic setup assuming the conduction and 
screening in the leads to be more efficient than in the 
conductor) . In addition, we assume that conduction elec- 
trons are completely polarized inside the ferromagnets, 
i.e. the population of carriers with a spin direction op- 
posite to that of a magnet is fully depleted (half-metallic 
ferromagnets). Therefore, when the polarizations of the 
leads are antiparallel, a conduction electron cannot be 
transferred across the valve without changing its spin 
direction. As a result the resistance of a spin-filter is 
very large unless there is a substantial amount of spin- 
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flip scattering inside the N-region. We assume the F-N 
interfaces to be spin-conserving but allow for the finite 
contact (tunnel) resistances R. 

Stochastic diffusion equations. The electron motion 
inside the N-region is diffusive with the mean free path 
much smaller than the size of the valve L (but yet much 
larger than the Fermi wavelength). At temperatures low 
enough the inelastic (electron-phonon, electron-electron) 
scattering is suppressed (once the inelastic scattering 
length exceeds L). The electron distribution is there- 
fore almost isotropic in momentum space and can be de- 
scribed by the spin and energy-dependent distribution 
functions f a (x,e), with a = ± being a spin index: + 
corresponding to spin-up electrons and — to spin-down 
electrons. 

If the system is driven out of equilibrium (e.g. by apply- 
ing a voltage bias to the leads) , the distribution function 
becomes spatially inhomogeneous resulting in the electric 
current (we assume the cross-sectional area of the valve 
to be equal to unity), 

, , a df a {x, e) _ . . . . 

e ox 

where a = e 2 vD is the conductivity in the N-region, v 
is the density of states per single spin direction and D 
is the diffusion constant. The last term in Eq. is the 
stochastic Langevin source. It has zero expectation value 
and a correlator that similarly to the spinless case 0] is 
determined by the mean value of the electron distribution 
function, 

J a (x,e,t)Jp(x>,e>,ti) = 2o-5 a/3 Aj a (x, e) [1 -J a (x, e)] , 

where we have abbreviated A = S(x — x')6(e — e')6(t — t') 
and assumed no summation over the repeated indexes. 
The stochastic source J a is due to the random indepen- 
dent (i.e. Poissonian) events of spin-conserving scattering 
from disorder. 

The particle conservation implies a second relation be- 
tween the electric current and particle density (here- 
inafter we drop the arguments when it could not lead 
to confusion), 

The first term in the right-hand side accounts for the 
average particle flow between states with opposite spins 
due to spin-flip scattering (customary in treating spin- 
dependent diffusion problems 0). The spin-flip length 
L s is assumed to be much larger than the mean free path 
but no restrictions as to its relation to the size of the 
system L are imposed. The last term in Eq. (j5J is the 
Langevin source for the spin-flip scattering arising from 
randomness of a spin-flip process. It is similar to the 
stochastic terms for the fluctuations of the number of 



photons in disordered optical media [15(. Its second mo- 
ment is equal to the mean flow between states with dif- 
ferent spin directions, 

C(x, e, t)C(x>,e>, i 7 ) = ^ £ 7 Q (1 - ~ a ). (3) 

s a 

which utilizes the fact that spin-flip scattering events 
are independent and obey Poissonian statistics. In writ- 
ing Eqs. (12i;j[) we suggested that the spin-flip scattering 
is energy-conserving. This assumption is well justified 
whenever a typical energy change during a spin-flip is 
small compared to the characteristic scale of the electron 
distribution (set by the temperature T or external bias 
eV). 

The above equations must be supplemented with ap- 
propriate boundary conditions. We assume that the in- 
terface resistances at the left and right contacts are the 
same R. Since there is no charge accumulation in the 
system, the diffusive currents should match the tun- 
neling currents through the interfaces. In particular, for 
the antiparallel valve configuration the boundary condi- 
tions read, 

3- = 7nUl- /-] + Tl, 3+ = 0, at x = -f 

3+ = Mf- - fn]+ J n, 3- = 0, at x = § . 1 1 

For the parallel configuration one has to interchange + 
and — indexes in the second line of Eq. Q). The stochas- 
tic sources Tl and Xr accounting for the randomness of 
the electron tunneling through the interfaces have (at 
T = 0) the variance [18| . 

Xi{e,t)T k (e',t') = S ik S(e-e')S(t-t')eJ(e), (5) 

here J(e) = ^2 a j a {x,e) is the total current independent 
of the coordinate x, as readily seen from Eq. (J2J. The 
current at the contacts is due to electrons with a single 
spin direction only. 

It is convenient to use the particle density and 
spin density distributions as well as the corresponding 
Langevin sources, 

fjs = \{f+±f-), J,Js = \(J + ±J-)- 

Combining Eqs. and J3J) we obtain (in the stationary 
regime) the equations for the particle and spin distribu- 
tion, 

(P£ = edJ ^ = A + ^ + £^ ( 6 ) 
dx 2 a dx ' dx 2 L 2 D a dx 

Note, that different Langevin terms (I, J ', C) are inde- 
pendent and have zero cross-correlators. 

Average electric current. The mean (averaged over 
time) solution of the equations @ with the boundary 
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conditions is elementary and yields the distribution 
function, 

7» = fj ^- /«] ( J " (*)) . (7) 

with R s = L s /a standing for the characteristic resistance 
on a spin- flip length L s . The total resistance Ro and the 
function M(x) depend on the magnetization of the leads. 
For the antiparallel configuration, 

Ro = R N + 2R + R s coths, M(x) = cosh ( x J L ^ ; (g) 



sinh ; 



while for the parallel configuration, 



Ra = R N + 2R + R s tanhs. Mix) = hmh ( x / Ls \ / Q \ 

coshs 

Here s = L/2L S is the dimensionless measure of the 
amount of spin-flip scattering in the system, and Rn = 
L/2o~ is the resistance of the normal region. 

The total mean electric current calculated with the 
help of Eqs. and (J7J) is determined by the total re- 
sistance, 

-J = Y^jde-Ue) = ±-f d e[f L -f R] = ^ (10) 

where the bias eV is the difference in the chemical po- 
tentials of the left and right leads /^(e — eV) = /n(e). In 
the absence of spin-flip scattering s — > the resistance of 
the parallel (the valve switched 'on') configuration tends 
to the 2R + 2Rn value, while for the antiparallel one 
(the valve switched 'off') it diverges. For s s> 1, both 
resistances tend to 2R + Rn- 

Shot noise. To solve Eqs. (JHJ it is convenient to write 
the fluctuating part of the distribution function in the 
form, 

Sf(x) = A + Bx + - [ dx'G (x,x')^- (11) 
a J ox' 

Sf s (x) — A s cosh (x/L s ) + B s sinh (x/L s ) 

+ /^O.^, (12, 

with the help of the Green function vanishing at the in- 
terfaces, 

_ sinh (x < /L s + s) sinh (x > /L s - s) 
Gs(X ' Xj "^ s snflT(2^ ' [U) 



with x<(x>) standing for the smaller (larger) of the two 
coordinates x,x'. The function Gq(x,x') is determined 
from the same expression (II 3D with s — ► 0. The coeffi- 
cients A, A s , £?, £? s are to be determined from the bound- 
ary conditions Q. It should be pointed out that the dis- 
tributions in the leads do not fluctuate Sfi = S/r = 0. 
The fluctuation of the total current is determined by the 
coefficient B only, according to, 



2a 2 f 

SJ = B + - dx J(x). 

e L J 



(14) 



Resolving a set of linear algebraic equations (obtained 
from the boundary conditions) with respect to B we find 
the fluctuation of the total energy- resolved current, 

SJ(e,t) = ^-Pl+M + ^tY; f dx K a (x)J a (x) 
Ho K L J 

dx M{x)£(x), (15) 



eR a D 



where the kernel function K a {x) depends on the valve 
configuration, 



KJx) 



a 



sinh (x/Ls) 

sinh s ' 
cosh (x / L s ) 

cosh s 5 



antiparallel, 
parallel. 



(16) 



The static shot noise power determined as the zero- 
frequency transform of the current-current correlation 
function S = J dt (5J(t)5J(0)) can now be calculated 
from Eq. Ijl5(l with the help of the correlation functions 
for the Langevin sources, 



s = ej2 -Rj + R%T.j dedx K(*)f«Q /«) 

+M 2 (x)U(l-J^)}. (17) 



Substituting the mean distribution functions Q into Eq. 
(|17fl and evaluating the spatial integrals we obtain the 
final expressions for the dimensionless noise-to-current 
ratio, F = S/eJ, also known as the Fano factor, 
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F, 



coth s 



,r 



2p 2 



11 



2P1T 
f tanh s 



2p 3 
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s[5 - cosh (4s)] + 2 sinh (2s) s 



2p 2 



II 



2p 3 



8 s sinh s 
s[5 - cosh (4s)] - 2 sinh (2s) 
8s cosh 4 s 



— — s coth s 



— s tanh s 



(18) 
(19) 



with p = Rq/R s being the dimensionless total resistance: 
pjj = s(r + 1) + coths for the antiparallcl configuration 
and pn = s(r + 1) + tanhs for the parallel configura- 
tion. We also introduced the dimensionless tunneling re- 
sistance r = 2R/Rn- 

Figs. 2 and 3 illustrate the Fano factor behavior with 
the spin-flip intensity s for different values of the contact 
resistance r for antiparallel and parallel valve configura- 
tions respectively. Let us first discuss the regime of trans- 
parent F-N interfaces r = 0. For large spin- flip scatter- 
ing, s — > oo, the shot noise approaches the universal value 
F = 1/3 independent of the relative magnetization of the 
leads. This is obvious since an injected electron quickly 
loses its polarization. For intermediate values, s > 1, 
the noise is slightly increased by spin-flip scattering both 
for parallel and antiparallel spin valve configurations. For 
small spin- flip intensity, s < 1, the noise behavior is com- 
pletely different. In the parallel configuration the Fano 
factor is returned to its universal value 1/3, which is easy 
to understand by realizing that electric current is trans- 
ferred predominantly by the spin-down states. In the 
antiparallel configuration, however, the small-amount of 
spin-flip scattering is responsible for the finite conduc- 
tance itself. The spin-flip induced fluctuations contribute 
to the noise comparably to the disorder-induced fluctu- 
ations. The noise power is therefore enhanced reaching 
ultimately the full Poissonian value usually reflective of 
the independent electron transmission, like in a tunnel 
junction or a Shottky vacuum diode. 

The presence of contacts with the finite resistance r 
changes the noise-to-current ratio. For large spin-flip 
scattering, 
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F 



I 



I 1 



+ 



1 



2(r + l) 2(r + l) 2 6(r + l) 3 



the Fano factor is increased monotonously from F = 1/3 
to F = 1/2 by changing r from zero to infinity. Exactly 
opposite, however, happens for antiparallel configuration 
with low spin-flip scattering ('off '-state of the valve), s < 
1, where the presence of contacts actually suppresses the 
noise power. 

The stochastic diffusion equations presented here allow 
for the discussion of the time-dependent problems as well, 
e.g. frequency dependence of the noise power. Without 
spin-flip scattering the noise spectrum is white as a result 
of the Debye screening . Shot noise in a spin valve is 
different since fluctuations of spin density do not require 



fluctuations of charge density. Mathematically it is illus- 
trated by the existence of the new (spin-flip) frequency 
scale D/L 2 S . The calculations would be similar to those 
performed for the phononic noise spectrum [l6| . 
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This formula holds for T = only. To consider the ef- 
fects of thermal fluctuations one needs a more general 
expression fo r the intrinsic noise of a tunnel barrier, 
l(e, t)I{e',t') = 6{e-e')5(t - 1')^ 1 [n L (l -n R ) + n H (l- 
fir,)], where ul and ur are the distribution functions at 
the left and right side of the barrier. 
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FIG. 2: Fano factor versus the spin-flip intensity for 
different values of the contact resistance: r — (solid line), 
r = 1 (dashed), r — 3 (dotted), r = 10 (dot-dashed). 



FIG. 3: Fano factor F±± versus the spin-flip intensity for 
different values of the contact resistance: r = (solid line), 
r = 1 (dashed), r — 3 (dotted). 



